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We present a measurement of the ratio of top quark branching fractions 
R = B(t — > Wb)/B(t — > Wq), where q can be a d, s or b quark, in the lepton+jets and dilepton 
tt final states. The measurement uses data from 5.4 fb _1 of pp collisions collected with the DO 
detector at the Fermilab Tevatron Collider. We measure R = 0.90 ± 0.04, and we extract the CKM 
matrix element |Vtf,| as |V t b| = 0.95 ± 0.02, assuming unitarity of the 3x3 CKM matrix. 

PACS numbers: 12.15.Hh, 13.85.Qk, 14.65.Ha 



The standard model (SM) of particle physics contains 
three generations of quarks. The top quark belongs to the 
third generation, and is of interest not only because of its 
large mass [l|, but also because its decay has not been 
examined in great detail, and may prove to be inconsis- 
tent with the SM. The decay rate of the top quark into 
a W boson and a down- type quark q (q = d, s, b) is pro- 
portional to |Vt q | 2 , the squared element of the Cabibbo 
Kobayashi Maskawa (CKM) matrix 2]. Under the as- 
sumption of a unitary 3x3 CKM matrix, \Vtb\ is highly 
constrained to \V tb \ = 0.9991521°;°°°°^° [|, and the top 
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quark decays almost exclusively to Wb. The existence 
of a fourth generation of quarks would remove this con- 
straint and accomodate significantly smaller values of 
\Vtb\- A smaller value of \Vtb\ could be observed directly 
through the electroweak production of single top quarks, 
for which the cross section is proportional to \V tb \ 2 , and 
could also affect the decay rates in the tt production 
channel. The latter can be used to extract the ratio of 
branching fractions R: 

_ B(t — s- Wb) = I Vtb | 2 m 

B(t -> Wq) I V a | 2 + \V ts \ 2 + \ V td \ 2 ' [ ' 

Given the constraints on the unitary 3x3 CKM matrix 
elements, R is expected to be 0.99830+°;°°°°^. Along with 
a measurement of | Vtb | using single top quark production, 
the measurement of R provides the possibility of a study 
of I- 

This Letter presents a measurement of R using a 
data sample corresponding to an integrated luminosity of 
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5.4 fb _1 of pp collisions, collected with the DO detector at 
the Fermilab Tevatronpp Collider at y/s = 1.96 TeV. We 
present measurements in the lepton+jets (£+]ets) chan- 
nel, in which one W boson from it — > W + qW~q produc- 
tion decays into a quark and an antiquark and the other 
into a charged lepton and a neutrino, and in dilcpton (ti) 
final states, in which both W bosons decay into Iv. We 
also present the combination of these two measurements. 
We consider events in which the charged leptons are ei- 
ther electrons or muons, produced directly from the W 
decay or from the leptonic decay of a r lepton. The re- 
sult from the £+jets channel corresponds to an improve- 
ment of the measurement using 0.9 fb _1 in which we 
extracted R > 0.79 at a 95% CL. This is the first DO 
measurement of R in the 11 channel. The CDF Collabo- 
ration has measured R in the ^+jets and tt channels in 
160 pb _1 of integrated luminosity [f|, and found a limit 
of R > 0.61 at 95% CL. 

Our measurement is performed by distinguishing be- 
tween the standard decay mode of the top quark ti — ¥ 
W + bW~b (indicated by bb), and decay modes that in- 
clude light quarks (qt = d,s): tt — > W + bW~q~i (bqi) and 
ti — > W + qiW~q~i (qiqi). The selection of an enriched tt 
sample and identification of jets from b quarks are crucial 
elements of the analysis. 

The DO detector [7J has a central tracking system con- 
sisting of a silicon microstrip tracker and a fiber tracker, 
both located within a 1.9 T superconducting solenoidal 
magnet, desi gned to optimize tracking at pseudorapidi- 



ties 1 77 1 < 3 21[. The liquid- argon/uranium calorime- 



ter has a central section covering pseudorapidities \rj\ up 
to s» 1.1, and two end calorimeters that extend cover- 
age to |t7| ps 4.2 8]. The outer muon system, covering 
|?7| < 2, consists of a layer of tracking detectors and scin- 
tillation trigger counters in front of 1.8 T toroids, fol- 
lowed by two similar layers behind the toroids Q . In the 
^+jets channel, we rely on the event selections used for 
the measurement of the tt production cross section [Icj . 
Details on object identification and selections are only 
briefly summarized as follows. We select ti events by 
taking advantage of their distinct topology. We require 
at least three jets within |7j| < 2.5, with transverse mo- 
mentum pt > 20 GeV, of which at least one has to 
have pt > 40 GeV. We require one electron (muon) 
of pt > 20 GeV, \rj\ < 1.1 (|jy| < 2.0) isolated from 
jets. In addition, events with a second isolated elec- 
tron or muon of pt > 15 GeV are removed in order to 
ensure that the t+jets and it samples are statistically 
independent. The imbalance in transverse energy, $t, 
must fulfill $ T > 20 GeV (fir > 25 GeV) in the e+jets 
(/i+jets) channel. The most important background in 
the i+'jets channel is from VF+jcts events which can pro- 
duce a similar final state to tt events. There is also sig- 
nificant background contribution from multijet produc- 
tion, in which a jet is misidentified as an electron, or a 
muon from the semileptonic decay of a hadron appears 



isolated. Smaller background contributions arise from 
electroweak single top quark production, Drell-Yan and 
Z boson production (decaying to Z + Z~+jets) or diboson 
production (WW, WZ or ZZ). The multijet background 
is estimated from control samples in data [Io| . while the 
tt signal and electroweak backgrounds are simulated us- 
ing Monte Carlo (MC) event generators ALPGEN and 
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12j . and a GEANT-based [13| simulation of 
the DO detector. Drell-Yan and Z boson production is 
normalized to the next-to-next-to-leading order (NNLO) 
QCD prediction 14j. All other electroweak backgrounds 
are normalized to their next-to-leading order (NLO) cross 
sections, while the M^+jets background is normalized to 
data using an iterative procedure (10|. 

For the 11 channel, we use the same selections as used 
for the measurement of the ti cross section described in 
Ref. 15j. In this final state, the ti signature consists of 
two energetic, oppositely charged isolated leptons, large 
JJst and two high pt jets. We consider separately the 
three final states ee, /i/i and e\x. For the e/i final state, we 
also consider events with only one reconstructed jet. To 
select ti events, we require two isolated leptons, electrons 
or muons, with pt > 15 GeV, |»y| < 1.1 or 1.5 < \rj\ < 2.5 
(| 77 1 < 2) for the electron (muon), and at least two jets 
with pt > 20 GeV and \rj\ < 2.5. For fifi events we 
require $t > 40 GeV. In the e^ channel, the sum of the 
transverse momenta of the lepton and two jets of highest 
Pt must be larger than 110 GeV. That sum must be 
higher than 105 GeV when only one jet is reconstructed. 
In the ee and fXfj, channels we use the $t significance to 
differentiate events with true $t from escaping neutrinos 
and events with $t arising from mismeasurement. The 
|?t significance for each event is defined in terms of a 
likelihood discriminant constructed from the ratio of $t 
to its uncertainty [16]. The significance is required to 
correspond to more than five. The main background in 
the ££ final states is composed of Drell-Yan, Z boson 
production and diboson events, and is estimated using 
MC simulation, normalized to the NNLO and NLO cross 
sections respectively. There is also a background from 
multijet events that we estimate from data [15j. 

We use a neural network (NN) 6-tagging algorithm 
to identify jets that contain b quarks, and thereby distin- 
guish the bb, bqi and qiqi final states. The inputs to the 
NN include impact parameters of tracks associated with 
the jets, and the properties of secondary vertices within 
the jet. Only taggable jets, i.e., jets matched to a set of 
tracks, are considered by the NN. For each taggable jet, 
we obtain an output from the NN which ranges between 
zero and twelve, with larger values more likely to corre- 
spond to jets originating from b quarks. Non-taggable 
jets are assigned the NN output value —1. 

We pursue different strategies to measure R in the 
t+jets and l£ channels. In the i'+jets channel, we count 
the number of jets that pass our threshold on the b- 
tagging NN output; this requirement has an efficiency 
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for b jets of 55 ± 4%, while admitting 1.5 ± 0.1% of light 
jets. The events are divided into subsamples according 
to lepton flavor (e or //), the number of jets in the event 
(3 and > 3), the data taking period (the first 1 fb _1 and 
the rest [10|), and the number of identified b jets (0, 1 or 
> 1). The events are separated further using a multivari- 
ate kinematic discriminant in subsamples dominated by 
background, i.e., events with zero 6-tagged jets, or one b- 
tagged jet in the sample with exactly three jets, and zero 
6-tagged jets in the sample with more than three jets. 
This discriminant is based on a multivariate technique 
(random forest of decision trees [lH) and uses several 
variables that exploit the kinematic differences between 
ti signal and background. In addition to ti MC sam- 
ples with SM decay ti — > WbWb, samples for the decay 
modes including light quarks are generated with pythia 
(ti — > WbWqi and ti — > WqiWqi), for a top quark mass 
of m t — 172.5 GeV. The expected number of ti events 
with m ^-tagged jets can be written as: 

fJ%(R,att) = [R 2 e m {bb) + 2R(1 - R)e m {bqi) 

+ (1 - R) 2 e m { mi )] a tf B 2 (t -> Wq)L , (2) 

where e m is the product of the selection efficiency and 
the probability of an event to have m 6-tagged jets for 
each of the three (bb, bqi and qiqi) decay modes, a t i is 
the ti production cross section and L is the integrated 
luminosity. A maximum likelihood fit is performed using 
the function: 



"£+jets 



Y[ P[n\ fi l {R, <T t h Vk)]P[tf M Ji Au\ x 

J]^K;0,sd) (3) 



i=l 



where i runs over the subsamples and bins of the multi- 
variate discriminant, and P[n, fi(R,a t t,iSk)] is the Pois- 
son probability to observe n events for an expected num- 
ber of /i(-R, <Jtii v h) events. The expectation \i(R, Utu v k) 
is the sum of the expected number of ti — > bb, bqi 
and qiqi events and the expected number of background 
events. The observed and expected numbers of multi- 
jet events are denoted n l MJ and \i l M j, and the Poisson 
terms P[n\ij,n l MJ ] take into account the fluctuation of 
the number of multijet events within the statistical un- 
certainties with which it is determined in dedicated data 
samples. Figures Q] (a) and (b) show the number of b- 
tagged jets in £+jets events for data and simulation for 
R = 0, R = 0.5 and R = 1. To reduce the dependence 
of the measurement on the input ti cross section, we si- 
multaneously extract <j t i from data, taking into account 
the three channels ti — > bb, bqi and qiqi. A parameter 
that accounts for each independent source of systematic 
uncertainty k is modeled by a Gaussian function Q with a 
mean of zero and a width corresponding to one estimated 
standard deviation (SD) of that uncertainty. This proce- 
dure correlates systematic uncertainties among channels 



by using the same parameter for a common source of 
systematic uncertainty. 

In the dilepton channels ee, fifj,, and e/i with at least 
2 jets, we apply the NN 6-tagging algorithm to the two 
jets of highest-pT, and use the smaller of the two NN 
outputs to calculate the likelihood, as it yields the best 
expected precision on R for values close to unity. The b- 
tagging algorithm is applied to the single reconstructed 
jet in the e/i channel with exactly 1 jet. We construct the 
templates for the decay modes bb, bqi, qiqi for ti as well as 
for all background components, forming the likelihood by 
running the product of Eq.[4]over all fourteen bins of the 
NN discriminant in all four channels, yielding thereby a 
product with 56 factors: 

C U= Y[Pln\tS(R,Vti,Vk)}Y[G(vk;0,SD) . (4) 

i=l k 

The expected number of events, ^(R,cr t t,Vk), is given 
by Eq. [2 where e m describes now the efficiency for the 
discriminant bin m, and v% can affect the individual com- 
ponents of ^ (R, (j t t) ■ Figure Q] (c) compares the distri- 
butions of the discriminant for predicted and observed 
events in the combined 11 final state. 

Several systematic uncertainties can impact the mea- 
surement of R. We consider the same sources of system- 
atic uncertainties as for the cross-section measurements 
in the ^+jets and it channels, and refer to Refs. 2, 15 1 
for details. The main source of systematic uncertainty 
on R is from the 6-tagging probability. Other important 
contributions to the systematic uncertainty on R arise 
from the jet identification efficiency, jet energy scale and 
resolution, and uncertainties on the background normal- 
ization as well as on modeling of the signal. The latter 
includes contributions from higher order effects, color re- 
connection, choice of parton distribution functions and 
initial and final-state gluon radiation. For consistency 
with Refs. [13, [H| we also quote separately the smaller 
systematic contributions from limited number of events 
in the templates and the uncertainties on the heavy-flavor 
fraction for the W+jets process, the trigger efficiency and 
lepton identification. We account for the fact that un- 
certainties from jet identification, jet energy scale and 
resolution, 6-jet identification, and higher-order correc- 
tions can affect the distribution of the discriminants in 
the £+jets channel, and the NN discriminants in the ££ 
channel. We verify that the measurement of R does not 
depend on m t by generating MC samples at different mt 
values. In the £+jets channel we obtain: 
R = 0.95 ± 0.07 (stat+syst) 

u t i = 7.90±°;$ (stat+syst) pb, 
and in the 11 channel 

R = 0.86 ± 0.05 (stat+syst) 

ott = 8.19tlst (stat+syst) pb. 
The results are in agreement with each other, and the 
extracted cross sections are consistent with those from 
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FIG. 1: (color online) (a) Number of 6-tagged jets in £+jets events with three jets and (b) at least four jets, (c) Distribution 
in the minimum fa-tag NN output of the jets of highest-pr for dilepton final states. 



Refs. [Tfl, 

In these a t i measurements, we do not as- 
sume that Bit — > Wb) = 1 as was done for the results 
in Refs. [HJ, , but only require Bit -> Wq) = 1. The 
combined measurement is obtained by fitting simultane- 
ously all channels in the 11 and ^+jets final states. This 
yields: 

R = 0.90 ± 0.04 (stat+syst) 
a tl = 7.74+°;^ (stat+syst) pb. 
Table [J summarizes the systematic uncertainties for the 
three results on R. While in the M channel the statisti- 
cal uncertainty still dominates, the £+jets and the com- 
bined result are dominated by systematic uncertainties. 
If we assume unitarity of the CKM matrix, we extract 
\Vtb\ — 0.95 ± 0.02. Constraining the tt cross section to 
the SM value of 7.5±g;f pb yields R = 0.90 ± 0.04, 
identical within rounding errors to the result of the si- 
multaneous fit. 

Using the combined result, we extract intervals on R 
as well as on \Vtb\ from Eq. fT]). assuming unitarity of the 
3x3 CKM matrix. By applying the frequentist approach 
using the likelihood ratio ordering principle proposed by 
Feldman and Cousins (20j ] . we obtain the intervals in R 
as 0.82-0.98 and V tb as 0.90-0.99 at 95% CL. The ex- 
pected limits are R > 0.92 and V tb > 0.96 at 95% CL. 
Figuredshows the bands for 68%, 95% and 99.7% confi- 
dence limits on R. Our result is compatible with the SM 
expectation at the 1.6% level. At 99.7% CL, we obtain 
R > 0.77 and |Vf 6 | > 0.88. 

Without assumptions on the unitarity of 
the CKM matrix, we can write Eq. ([1]) as: 
(l-R)/R=(\V ts \ 2 + \V td \ 2 )/\V tb \ 2 , and set a limit 
on this ratio at 99.7% CL of: (1 -R)/R< 0.30. 

To summarize, we have measured the ratio of branch- 
ing fractions R = B(t — > Wb)/B{t -> Wq) in both lep- 
ton+jets and dilepton channels. In the combined analy- 
sis, we find R — 0.90 ±0.04, which agrees within approxi- 
mately 2.5 standard deviations with the SM prediction of 
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FIG. 2: (color online) Limit bands at 68%, 95% and 99.7% 
CL on R, with the measured value (dotted line). 



R close to one. This is the most precise determination of 
R to date. Using the approach of Ref. [20( and assuming 
the unitarity of the CKM matrix, we extract the interval 
at 95% CL on the element V tb as 0.90-0.99. 
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TABLE I: Uncertainties on the measurements of R in the U and £+jets channels as well as for the combination of the two. 
We evaluate the impact of each class of systematic uncertainties by calculating R and er t j using the corresponding parameters 
v shifted by ±1SD from their fitted mean. The final line shows the quadratic sum of the systematics, which can be slightly 
different from the one obtained with the global fit. 



it 


^+jets 


Combination 


Source 


+SD 


-SD 


+SD 


— SD 


+SD 


-SD 


Statistical 


0.041 


-0.042 


0.030 


-0.029 


0.023 


-0.023 


Muon identification 


0.002 


-0.002 


0.000 


-0.001 


0.001 


-0.001 


Electron identification and smearing 


0.004 


-0.004 


0.000 


-0.000 


0.001 


-0.002 


Signal modeling 


0.007 


-0.006 


0.009 


-0.011 


0.004 


-0.006 


Triggers 


0.003 


-0.003 


0.001 


-0.001 


0.002 


-0.002 


Jet energy scale 


0.008 


-0.008 


0.017 


-0.016 


0.003 


-0.008 


Jet reconstruction and identification 


0.010 


-0.009 


0.018 


-0.022 


0.009 


-0.013 


6-tagging 


0.018 


-0.019 


0.065 


-0.056 


0.034 


-0.033 


Background normalization 


0.020 


-0.020 


0.004 


-0.005 


0.008 


-0.010 


W fractions matching + higher order effects 






0.001 


-0.001 


0.001 


-0.002 


Instrumental background 


0.013 


-0.013 


0.003 


-0.004 


0.005 


-0.007 


Luminosity 


0.010 


-0.010 


0.001 


-0.001 


0.004 


-0.004 


Other 


0.002 


-0.002 


0.000 


-0.000 


0.001 


-0.001 


Template statistics for template fits 


0.002 


-0.002 


0.011 


-0.011 


0.010 


-0.010 


Quadratic sum of systematics 


0.035 


-0.035 


0.071 


-0.064 


0.038 


-0.040 
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